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, DAMA is an observatory for rare processes based on the development and use 

""^5 ' of various kinds of radiopure scintillators. Several low background set-ups have 

been realized with time passing and many rare processes have been investigated. 
Main activities are briefly summarized in the following and the main arguments 
on the the results achieved in the investigation of the WIMP annual modulation 
signature are addressed. Next perspectives are also mentioned. 
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ABSTRACT 



1. Introduction 

DAMA is an observatory for rare processes based on the development and use of 
various kinds of radiopure scintillators. The main experimental set-ups are: i) the 
~ 100 kg Nal(Tl) set-up, which has completed its data taking in July 2002; ii) the 
new 250 kg Nal(Tl) LIBRA (Large sodium Iodide Bulk for RAre processes) set-up, 
whose installation is started at fall 2002; iii) the ~ 6.5 kg liquid Xenon (LXe) pure 
scintillator; iv) the R&D installation for tests on prototypes and small scale experi- 
ments. Moreover, in the framework of devoted R&D for higher radiopure detectors 
and PMTs, sample measurements are regularly carried out by means of the low back- 
ground DAMA/Ge detector, installed deep underground since about a decade and, 
in some cases, at Ispra. 

The locations of the DAMA experimental installations in the Gran Sasso under- 
ground laboratory of I. N.F.N, are shown in Fig. ^ 

a also: Universita' "Campus Biomedico" di Roma, 00155, Rome, Italy 
b also: University of Zhao Qing, Guang Dong, China 
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Figure 1: The locations of the DAM A experimental installations in the Gran Sasso underground 
laboratory of I.N.F.N. 



2. DAMA/LXe 

The DAMA/LXe experiment has followed the former Xelidon experimen on R&D 
developments of liquid Xenon (LXe) detectors and has realized since '90 several LXe 
scintillator prototypes using natural Xenon. Then, it has preliminarily put in mea- 
cement the set-up used iu the data takiug of ret. HO by using Kr-free Xenon 
enriched in 129 Xe at 99.5%. This set-up was significantly upgraded at fall 1995 (as 
mentioned e.g. in ref. "*) by including: i) a new purification system without Ox- 
isorb; ii) a new low background Pb in the shield; iii) a new low background Cu in the 
shield and in the insulation vessel; iv) a substitution of some PMTs and improved low 
background voltage dividers, etc. As a consequence, a significant improvement in the 
counting rate (as tipically experienced in low background experiments; see e.g. the 
case of ref. ^ for Ge detector) was obtained. The rate measured after this upgrading 
was well consistent over the whole energy spectrum and during several years of data 
taking, although some other minor substitutions have obviously been carried out with 
time passing. 

In summer 2000 the set-up was again deeply modified (reaching the configuration 
reported in Fig.5 of r ef. j^ ) to handle also Kr-free Xenon enriched in 136 Xe at 68.8% 
and in 134 Xe at 17.1% Ed. The main differences among the previous and the present 
experimental set-up operating with this latter gas are: the gas c , part of the set-up d 

c It was produced by another factory and many years before than the one enriched in 129 Xe; in 
addition, it was also previously used in a different underground experiment where different materi- 
als and vacuum/purification/filling/recovery system were operative and, afterwards, stored under- 
ground for long time in bottles with different possible effects from surface degassing, etc.; thus, the 
background is significantly different. 

d The vacuum/purification/filling/recovery system has been significantly modified to allow the 



and the isotope 6 . In particular, the high energy rate measured with this new set-up 
has been found higher than the one previously measured with the 129 Xe set-up. 

The main features of the set-up (see Fig. EJ) are described in ref. ^ and in related 
papers. We also remind that careful neutron calibrations have been carried out during 




Figure 2: On the left: the vessel of the LXe set-up. On the right: the vacuum/purification/filling/re- 
covery system and the passive shield ahead. 

several years IBE3_ 

As regards the more recent results achieved with this set-up on the Dark Matter 
investigation we mention the limits on recoils obtained by investigating the WIMP- 
129 Xe elastic scattering exploiting the pulse shape discrimination technique ^ and 
thos e obtained - in a given model framework - on the WIMP- 129 Xe inelastic scattering 

The same experiment has allowed to investigate several other rare processes such 
as nuclear level excitation of 129 Xe during charge-non-conserving processes and 



the possible electron decay through the channel: e 
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processes 
. In addition, 



the 



nucleon and di-nucleon decay into invisible channels has been investigated with a new 
approach based on the search for the radioactive daughter nuclei, created after 
the nucleon or di-nucleon disappearance in the parent nuclei. The advantage of this 
approach is a branching ratio close to 1 and an efficiency - since the parent and the 
daughter nuclei are located in the detector itself - also close to 1. Competitive limits 
have been obtained. 

After the latest upgrading of the set-up double beta decay modes in JffXe and in 
134 Xe have been deeply investigated reaching competitive limits as well . 



allocation and handling of the Xenon enriched in 136 Xe and in some other parts such as the cold 
trap (a new concept one) and the shield. 

e It can also be interested by different physical processes than the 129 Xe previously used. 



The data taking is continuing. 



3. DAMA/R&D 

The set-up named "R&D" is used for tests on prototypes and small scale experi- 
ments. A view of the passive shield of this installation is given in Fig. EJ 




Figure 3: View of the open shield of the R&D installation. 

This set-up which has been deeply upgraded in 2000/2001 has been used for mea- 
surements on low background prototype scintillators and PMTs realized in various 
R&D efforts with industries. Moreover, it is regularly used to perform small scale ex- 
periments mainly investigating double beta decay modes in various isotopes. Among 
the obtained results we remind the search for: i) (3(3 decay modes in 136 Ce and in 
i42 Ce |U. j.) 2EC2u decay mode in 40 CaE3; iii) (3(3 decay modes in 46 Ca and in 
40 CaE3; iv) (3(3 decay modes in 106 CdG3; v) (3(3 and/? decay modes in 48 CaHHJ; vi) 
2EC2v in 136 Ce and in 138 Ce and a decay in 142 Ce^. 

Fig. |U summarizes the results obtained in the searches for double beta decay 
modes. 
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Figure 4: Summary of the limits obtained by DAMA on various double beta decay processes. See 
text. 



4. DAMA/Ge for sampling measurements 

Various R&D developments to improve low background set-ups and scintillators as 
well as new developments for higher radiopure PMTs are regularly carried out. The 
related measurements on samples are performed by means of the DAMA low back- 
ground Ge detector f , which is operative deep underground in the low background 
facility of the Gran Sasso National Laboratory since many years. 



5. DAMA/Nal 



The main goal of the ~ 100 kg Nal(Tl) (DAMA/Nal) set-up has been the search 
for WIMPs by the annual modulation signature, which it has investigated over seven 
annual cycles; the data o f 4 of them (about 60000 kg-d) have been already released 
i z,u | z,i | z,z, | z,j | z,^ | z,^ | z,u | z, | z,c> ) anc j ma j n ar g umen ts will be discussed in the following. 

For completeness, we remind that - in addition to the investigation of the WIMP 
component in the galactic halo by means of the annual modulation signature - other 
approaches have also been exploited with DAMA/Nal such as the pulse shape discrim- 
ination technique and the investigation of possible diurnal effects Moreover, 
also ex otic Da rk Matter candidates such as neutral SIMPs, neutral nuclearities and 
Q-balls l ^l l ^ l have been searched for. A devoted search for solar axions has been 
carried out as well (see Fig. EJ) 
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Figure 5: Exclusion plot in the plane axion to photon coupling constant, g ail , versus axion mass, 
m a , achieved by DAMA/Nal in ref. 33) . The limit quoted in the paper (gt a77 < 1.7 x lO^GeV^ 1 
at 90% C.L.) is shown together with the expectations of the KSVZ and DFSZ models; see ref. 33 ) 
for details. 

In addition, DAMA/Nal has also allowed to investigate several other rare processes 
such as e.g. Pauli exclusion principle violation with spontaneous emission of protons 

This detector was specially realized with a low Z window to be sensitive to external radiation 
down to about ten keV. 



in Na and I nuclear level excitation of I and Na during charge-non- 
conserving processes and electron stability and non-paulian transitions in Iodine 
atoms (by L-shell) ^EJ. 

The set-up and its performances have been described in details in ref. ^^; since then 
some upgrading has been carried out. In particular, in summer 2000 the electronic 
chain and data acquisition system have been completely substituted, while during 
August 2001 the new HV power supply system and the new preamplifiers prepared 
for the foregoing LIBRA set-up have been installed here. 

The set-up has completed its data taking in July 2002. Various kind of data analyses 
are continuing; in particular, the total statistics of 107731 kg-d will be presented in 
near future giving further results on the WIMP annual modulation signature. 



5.1. Results on the WIMP annual modulation signature 

DAMA/Nal has been proposed in '90 and has been realized to investigate 
with suitable sensitivity the WIMP annual modulation signature, originally proposed 
in ref. ES. It has been so far the only experiment able to test this model indepen- 
dent signature for WIMPs with suitable mass, sensitivity and control of the running 
parameters. 

This WIMP model independent signature is based on the annual modulation of 
the signal rate induced by the Earth revolution around the Sun; as a consequence, 
the Earth will be crossed by a larger WIMP flux roughly in June (when its rotational 
velocity is summed to the one of the solar system with respect to the Galaxy) and 
by a smaller one roughly in December (when the two velocities are subtracted). The 
fractional difference between the maximum and the minimum of the rate is of order 
of ~ 7%. Therefore, to point out the modulated component of the signal, large mass 
apparata with suitable performances and control of the operating conditions - such 
as the ~ 100 kg highly radiopure Nal(Tl) DAMA set-up - are necessary. 

The annual modulation signature is very distinctive since a WIMP-induced sea- 
sonal effect must simultaneously satisfy all the following requirements: the rate must 
contain a component modulated according to a cosine function (1) with one year pe- 
riod (2) and a phase that peaks around ~ 2 nd June (3); this modulation must only 
be found in a well-defined low energy range, where WIMP induced recoils can be 
present (4); it must apply to those events in which just one detector of many actually 
"fires", since the WIMP multi-scattering probability is negligible (5); the modulation 
amplitude in the region of maximal sensitivity must be <7% (6). Only systematic 
effects able to fulfil these 6 requirements could fake this signature; no one has been 



found nor suggested by anyone after several years of investigations. 



5.1.1. The model independent evidence by DAMA/Nal 



A model independent analysis of the data of the four annual cycles has offered an 
immediate evidence of the presence of an annual modulation of the rate of the single 
hit events in the lowest energy interval (2-6 keV) as shown in Fig. |3 The x 2 test on 
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Figure 6: Model independent residual rate for single hit events, in the 2-6 keV cumulative energy 
interval, as a function of the time elapsed since January 1-st of the first year of data taking. The 
expected behaviour of a WIMP signal is a cosine function with minimum roughly at the dashed 
vertical lines and with maximum roughly at the dotted ones. 



the data of Fig. |H1 disfavors the hypothesis of unmodulated behaviour (probability: 
4 • 1CT 4 ), while fitting these residuals with the function A- cosu(t — t ), one gets: i) 
for the period T = 2e = (1.00 ± 0.01) year when t is fixed at the 152.5 th day of the 
year (corresponding to ~ 2 June); ii) for the phase to = (144 ± 13) days, when T 
is fixed at 1 year. In the two cases A is: (0.022 ± 0.005) cpd/kg/keV and (0.023 ± 
0.005) cpd/kg/keV, respectively. Similar results, but with slightly larger errors, are 
found in case all the parameters are kept free. 

The modulation is absent in other part of the energy spectrum as discussed e.g. 
in ref. ^J; for example, in the 6-10 keV energy region just above a modulation 
amplitude A = -(0.0017 ± 0.0037) cpd/kg/keV is obtained. 

We have extensively discussed the results of the investigations of all sources of 
possible systematics when releasing the data of each annual cycle; moreover, a dedi- 
cated paperESJ h as been devoted to this subject. No systematic effect or side reaction 
able to mimic a WIMP induced effect (that is to be both quantitative relevant and 
able to satisfy the six peculiarities of the WIMP signature) has been found. This 
deep discussion will not be repeated here since the reader can directly find it in ref. 
^ with all the details. 

In conclusion, a WIMP contribution to the measured rate has been candidate 
by the result of the model independent approach independently on the nature and 



coupling with ordinary matter of the involved WIMP particle. This is the main 
experimental result of DAMA/Nal. 

No other experiment has been realized so far able to exploit this model inde- 
pendent approach with similar exposed mass, sensitivity and control of the running 
parameters as DAMA/Nal. 



5.1.2. The corollary model dependent analyses searching for the nature of a candidate 

A corollary investigation can be pursued to investigate the nature and coupling 
with ordinary matter of a possible candidate. In this suitable energy and time 

correlation analysis is necessary as well as the choice of a complete model framework. 
We remark that a model framework is identified not only by the general astrophysical, 
nuclear and particle physics assumptions, but also by the set of values used for all 
the experimental and theoretical parameters needed in the calculations (for example 
WIMP local velocity, Vq, form factors and related parameters, quenching factors, halo 
model, coupling, etc., which are affected by relevant uncertainties) g . 

For simplicity, initially we have considered the particular case of purely spin- 
independent (SI) coupled WIMP. In fact, often the spin-independent interaction with 
ordinary matter is assumed to be dominant since e.g. most of the used target-nuclei 
are practically not sensitive to SD interactions as on the contrary 23 Na and 127 I 
are and the theoretical calculations are even more complex when including also this 
latter kind of intera ction. Moreover, only one set of possible parameters values has 
been adopted 1 ^0 1 22 1 Then, this case has been extended by c onside ring some of the 



existing uncertainties on the astrophysical velocity distribution 1 g i t and the physical 
constraint which arises from the upper limits on recoils measured by the same set-up 
More recently an investigation on the effect induced on the result by considering 
other possible and consistent halo models still for the particular case of purely SI 
coupled WIMPs has also been carried out in ref. 

Moreover, some of the other possible scenarios have also been considered such 
as extensions to the general case of WIMPs with both spin-independent (SI) and 
spin-dependent (SD) coupling and to the case of WIMPs with preferred inelastic 
scattering 1^3. In these latter cases, the effect of the uncertainties on some of the 
parameters has been included, but still only the simplified, approximate and non- 
consistent isothermal halo model has been considered in the calculations. 

Theoretical implications in terms of a neutralino with dominant SI interaction 
have been discussed e.g. in ref. I ^2 | 43 t f or some theoretical model frameworks and in 
terms of an heavy neutrino of the fourth family in ref. 



s Note that results given in terms of exclusion plots by experiments such as BSEOBl a i so d e _ 
pend on their own choice of nuclear, particle and astrophysical assumptions and of experimen- 
tal/theoretical parameters values; thus they have no generality at all. 



5.1.2.1 WIMPs with dominant SI interaction in given model frameworks 



A full energy and time correlation analysis - properly accou nting f or the phys- 
ical constraint arising from the measured upper limit on recoils l ^5 | 29| _ nag b een 
carried out in the framework of given model for purely spin-independent coupled can- 
didates with mass above 30 GeV (this bound being inspired by the lower bound on 
the supersymmetric candidate, as derived from the LEP data in the usually adopted 
super symmetric schemes based on GUT unification assumptions). A standard max- 
imum likelihood method has been used. Note that different model frameworks (see 
above) vary the theoretical expectations and, therefore, the best fit values of cross 
section and mass (as well as the allowed region) also vary. In particular, the inclusion 
of the uncertainties associated to the models and to every parameter in the models 
themselves as well as other possible scenarios largely enlarges the allowed region as 
discussed e.g. in ref. E31 f or the particular case of the astrophysical velocities and 
offers very large sets of best fit values. We take this occasion to stress that ref. EH 
also accounts for several models; thus, for several sets of best fit values h . Therefore, 
claims for contradiction in model dependent comparisons by "choosing" a particular 
set of best fit values, given there as an example, are - also in this respect - arbitrary 
and substantially wrong. 

As mentioned, more recently possible departures from the isothermal sphere model, 
which is the parameterisation usually adopted only in direct WIMP searches to de- 
scribe the halo although approximate and non-consistent, have been investigated in 
a systematic way in ref. ^* by some of us in collaboration with N. Fornengo and 
S. Scopel. Modifications arising from various matter density profiles, effects due to 
anisotropies of the velocity dispersion tensor and rotation of the galactic halo have 
been specifically investigated. In particular, several halo models with potential and 
matter density having a spherical symmetry have been investigated as well as halo 
models with spherical symmetry but anisotropic WIMP velocity distribution and halo 
models with axial symmetry (in these latter cases possible co-rotation or counter- 
rotation of the dark halo has also been considered). Also some triaxial models have 
been investigated. 

The global results of this analysis - which includes also those presented in ref. El 
- are shown in Fig. where the allowed regions obtained for the various considered 
halo models in given framework for three of the possible values of the local velocity, 
Vq, above WIMP mass of 30 GeV are superimposed. Obviously different best fit values 
correspond to each one. The cumulative result, which gives a direct impact of the 



h For example, for the particular model frameworks and assumptions of ref. where the WIMP 
local velocity, vq, has been varied from 170 km/s to 270 km/s to account for its present uncertainty, 
we obtained the best fit values mw = (72±i|) GeV and ^a S i = (5.7 ± 1.1) • 10~ 6 pb for v = 170 
km/s and m w = (43±g 2 ) GeV and £cr 5/ = (5.4 ± 1.0) ■ 10" 6 pb for v = 220 km/s. Here, £ is the 
WIMP local density in 0.3 GeV cm -3 unit, asi is the point-like SI WIMP-nucleon generalized cross 
section and mw is the WIMP mass. 
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Figure 7: Superposition of the regions allowed at 3<r C.L. in the given model frameworks by con- 
sidering the velocity distribution of each one of the halo models in ref. 28 ^. Three of the possible 
values of vo are considered. Obviously a specific set of best fit values for the WIMP mass and cross 
section corresponds to each region. Note that inclusion of other existing uncertainties will further 
enlarge the regions and increase the sets of best fit values. 

effect induced on the region allowed in the considered scenario only by the present poor 
knowledge of the right halo model, is shown in Fig. |HJ This region is compared with 
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Figure 8: Region allowed at 3a C.L. given by the superposition of all the allowed regions obtained, 
in the given model frameworks, considering several possible non-rotating halo models 2S \ Note that 
in these calculations only uncertainties on the halo model have been considered; the inclusion of the 
other existing uncertainties would further enlarge it. It is evident that this cumulative region, as 
those given e.g. in ref. 23 > 24 ) 5 accounts for a large set of best fit values for WIMP mass and cross 
section. The shaded region (which corresponds to the particular case of the approximate and non- 
consistent isothermal sphere halo model when assuming also vq = 220 km/s and po — 0.3 GeV/cm 3 
for the astrophysical, nuclear and particle physics assumptions and fixed parameters of ref. 28 ') is 
shown to point out only the effect due to the poor knowledge of the right halo model. 



the one obtained when considering, for a particular model framework (of the many 
possible), the approximate and non-consistent isothermal sphere halo model assuming 



in particular also vq = 220 km/s e po = 0.3 GeV/cm 3 . As one can see, the cumulative 
3 a C.L. allowed region is extended up to my/ ~ 270 GeV with cross section on 
nucleon in the range 10 _10 nbarn < £,o~si < 6 • 10~ 8 nbarn. Maximal co-rotating and 
counter-rotating models can extend the allowed region up to mw — 500 — 900 GeV 
(see Fig. 0). We further stress that in this analysis no other uncertainty than the halo 
model 1 has been considered; the proper inclusion of the other existing uncertainties 
will further extend the cumulative allowed region and offer further sets of best fit 
values. In particular, as already stressed also e.g. in ref. ^H, in all the analyses 
cautiously the Helm SI form factor has been adopted, which is the most cautious one 
for Iodine. The use of other form factors enhancing the expected signal from Iodine 
recoils would yield allowed regions corresponding to lower cross section values. The 
same would be if a spin-dependent component different from zero would be introduced 



5.1.2.2 WIMPs with mixed coupling in given model frameworks 

Since the 23 Na and 127 I nuclei are sensitive to both SI and SD couplings - on the 
contrary e.g. of nat Ge and nat Si which are sensitive mainly to WIMPs with SI coupling 
(only 7.8 % is non-zero spin isotope in nat Ge and only 4.7% of 29 Si in na 'Si) - the 
analysis of the data has been extended considering the more general case of a WIMP 
having not only a spin-independent, but also a spin-dependent coupling different from 
zero J . 

In this case the free parameters are £057, (,o- SD and m w for each given 6 value, 
where o~sd is the point-like SD WIMP cross section on nucleon and tgQ is the ratio 
between the effective SD coupling constants on neutrons, a n , and on proton, a p ; 
therefore, 9 can assume all the values between and n depending on the type of SD 
coupling. 

Note that the results in ref. I^SJ have been obtained by considering there only 
the approximate non-consistent isothermal sphere model to describe the galactic halo 
and a maxwellian WIMP velocity distribution with inclusion of some uncertainties 
on vq, on the nuclear radius and on the nuclear surface thickness parameter in the 
used Helm SI form factor, on the b parameter in the used used Ressel e al. SD 

'We note that, although the large number of halo models considered in this analysis, many other 
halo models are still possible, available and not yet considered here. 

•■We take this occasion to note that on the contrary of what is reported onE^, ref. ESJ does not 
exclude at all a possible spin-dependent solution. In fact, it refers only to the two particular cases 
(of the many possible) for purely spin-dependent coupling: specifically the case where the effective 
coupling constant on proton is equal to zero and the case where the effective coupling constant 
on neutron is equal to zero. In addition, not only these two cases are handled in strongly model 
dependent mode, but are even based on the wrong use of modulation amplitudes calculated for a 
par ticular purely spin-independent case. Thus no restriction at all arises - by the fact - from ref. 
Issil for any possible purely SD solution. In addition, the mixed case was not involved at all in that 
discussion. 




form factor and on the quenching factors. As an additional source of uncertainty we 
mention that an universal formulation is not possible for the SD form factor, thus 
other formulations than the one adopted here are possible and can be considered with 
evident implications on the obtained results. The same is for the used spin factor. 

For simplicity, Fig. |3] shows slices for some my/ of the region allowed - in the 
model frameworks of ref. - at 3 cr C.L. in the (£05/, C,ctsd, m w) space for four 
particular couplings: i) 9 = (a n =0 and a p ^ or \a p \ » \a n \); ii) 9 = it /A 
(a p = a n ); iii) 9 = n/2 (a n ^ and a p = or \a n \ » \a p \); iv) 9 = 2.435 rad = 
-0.85, pure Z° coupling). The case a p = —a n is nearly similar to the case iv). 
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Figure 9: A mixed SI/SD case: example of slices of the region allowed at 3 cr C.L. in the (£o\s/, 
$,<Jsd, mw) space for some mw and 9 values in the model frameworks considered in ref. 2& \ Only 
four particular couplings are reported here for simplicity: i) 9 = 0; ii) 9 = tc/4 iii) 9 = n/2: iv) 9 = 
2.435 rad. Note that e.g. Ge experiments are sensitive mainly only to SI coupling and, therefore, 
cannot explore most of the DAMA allowed regions in this scenario; the same is in most cases also 
for nat Xe since the odd-spin isotopes have the neutron as unpaired nucleon. These allowed regions 
would be further enlarged by taking into account the uncertainties existing on the halo models and 
related parameters, on the SD form factor (which has a not an universal formulation) and on some 
other experimental and theoretical parameters and assumptions. 



As already pointed out, when the SD contribution goes to zero (y axis in Fig. 
EJ), an interval not compatible with zero is obtained for £057. Similarly, when the SI 
contribution goes to zero (x axis in Fig. EJ), finite values for the SD cross section are 
obtained. Large regions are allowed for mixed configurations in the considered model 
frameworks also for £057 < 10~ 5 pb and £,<Jsd < 1 pb; only in the particular case of 
9 = I (that is a p = and a n ^ 0) (,<tsd can increase up to ~ 10 pb, since the 23 Na 



and 127 I nuclei have the proton as odd nucleon. Moreover, in ref. ^* we have also 
pointed out that: i) finite values can be allowed for £,<7sd even when £agi ~ 3 ■ 10~ 6 
pb as in the region allowed in the pure SI scenarios considered above; ii) regions not 
compatible with zero in the £<Jsd versus mw plane are allowed even when ^asi values 
much lower than those allowed in the dominant SI scenarios previously summarized 
are considered; iii) best fit values with both £057 and £,(Jsd different from zero are 
present for some mw and 9 pairs; the related confidence level ranges between ~ 3 a 
and ~ 4 cx^. 

Further investigations are in progress on these model dependent analyses to ac- 
count for other known parameters uncertainties and for possible different model as- 
sumptions. In fact, as mentioned, when including the uncertainties on the halo models 
and their parameters, on the SD form factor, on the spin factor and on some other ex- 
perimental and theoretical parameters, the allowed volume in the space (£crsi, ^sd, 
mw) for each 9 value would be further enlarged as well as the number of obtained 
sets of best fit values for the cross sections and the WIMP mass. 

In conclusion, this analysis has shown that the DAMA data of the four annual cy- 
cles, analysed in terms of WIMP annual modulation signature, can also be compatible 
with a mixed scenario where both £o"sj and £<tsd are different from zero. 

5.1.2.3 Inelastic Dark matter 

It has been suggested in ref. that the observed annual modulation effect could 
be induced by possible inelastic Dark Matter: relic particles that prefer to scatter 
inelastically off of nuclei. The inelastic Dark Matter could arise from a massive 
complex scalar split into two approximately degenerate real scalars or from a Dirac 
fermion split into two approximately degenerate Majorana fermions, namely x+ and 
X-, with a 5 mass splitting. In particular, a specific model featuring a real component 
of the sneutrino, in which the mass splitting naturally arises, has been given in ref. 
It has been shown that for the X- inelastic scattering on target nuclei a kinematical 
constraint exists which favours heavy nuclei (such as 127 I) with respect to lighter ones 
(such as e.g. nat Ge) as target-detectors media. In fact, X- can only inelastically 
scatter by transitioning to x+ (slightly heavier state than x~) and this process can 
occur only if the X- velocity is larger than v t hr = \J m^ N where mwN is the WIMP- 
nucleus reduced mass (c = 1). This kinematical constraint becomes increasingly 
severe as the nucleus mass, mjv, is decreased Moreover, this model scenario 

gives rise - with respect to the case of WIMP elastically scattering - to an enhanced 
modulated component, S m , with respect to the unmodulated one, Sq, and to largely 
different behaviours with energy for So and S m (both show a higher mean value) 

A dedicated energy and time correlation analysis of the DAMA annual modulation 
data has been carried out handling aspects other than the interaction type as in 
ref. EEl. In this scenario of Dark Matter with inelastic scattering an allowed volume 
in the space (£cr p , S, mw) is obtained!^. For simplicity, Fig. ITU1 shows slices of 
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Figure 10: An inelastic case: slices at fixed WIMP masses of the volume allowed at 3 a C.L. in 
the space (£cr p , 5, mw) obtained for the model frameworks considered in ref. 27 ); some of the 
uncertainties on used parameters have been included. Note that e.g. Ge experiments cannot explore 
most of the DAMA allowed regions in this scenario. These allowed regions would be further enlarged 
by taking into account the uncertainties existing on the halo models and their parameters and on 
some other experimental and theoretical parameters. 



such an allowed volume at some given masses (3 o C.L.) for the model frameworks 
considered in ref. . It can be noted that when 3> m^, the expected differential 
energy spectrum is trivially dependent on mw and in particular it is proportional to 
the ratio between ^a p and my/'-, this particular case is summarized in the last plot 
of Fig. The allowed regions have been obtained - as in the previous cases - 

by the superposition of those obtained when varying the values of some of the used 
parameters according to ref. Of course, each set of parameters' values gives rise 
to a different expectation, thus to different "most likely" values. Just as an example 
we mention that when fixing the other parameters as in ref. y the "most likely" 
values for a WIMP mass of 70 GeV are: i) £a p = 2.5 x KT 2 pb and 5 = 115 keV 
when Wo = 170 km/s, ii) £<t p = 6.3 x 10~ 4 pb and 5 = 122 keV when t> = 220 km/s; 
they are in 5 region where e.g. Ge and Si experiments are disfavoured. 

Finally, we note again that the allowed regions are further enlarged when properly 
including the uncertainties on the halo models, on the experimental and theoretical 



parameters and on the other assumptions. 



5.1.3. Comparison with some model dependent results 

5.1.3.1 ... from direct searches 

As mentioned above no other experiment directly comparable with the model 
independent DAMA/Nal result on WIMPs in the galactic halo is available at present; 
thus, claims for contradict ion are in trinsically arbitrary/wrong. 

Only few experiments EHEDEU^ which use different target nuclei and different 
methodological approaches, have released extremely poor selected statistics quoting 
an exclusion plot in a given particular model framework. 

Table ^ shortly summarizes some main items. As it can be seen, the mentioned 
exclusion plots are based on a huge data selection (releasing typically extremely poor 
exposures with respect to generally long data taking and, in some cases, to several 
used detectors). Moreover, their counting rate is very high and few/zero events are 
claimed after applying several strong and hardly safe rejection procedures (involv- 
ing several orders of magnitude; see Table HJ). These rejection procedures are also 
poorly described and, often, not completely quantified. Moreover, most efficiencies 
and physical quantities entering in the interpretation of the claimed selected events 
have never been discussed in the needed details. As an example, we mention the 
case of the quenching factor of the recoil target nuclei in the whole bulk material for 
the bolometer cases, which is arbitrarily assumed to be 1 k , implying a substantially 
arbitrarily assumed energy scale and energy threshold (and consequently arbitrary 
exclusion plot). Further uncertainties are present when, as in ref. a neutron 

background modeling and subtraction is pursued in addit ion. 

As regards in particular the Zeplin-I result of ref. EEH^ a very low energy 
threshold is claimed (2 keV), although the light response is very poor: between ~ 1 
ph.e./keV^^ (for most of the time) and ~ 2.5 ph.e./keV (claimed for 16 days) 
1 . Moreover, a strong data filtering is applied to the high level of measured counting 
rate (~ 100 cpd/kg/keV at low energy, which is nearly two orders of magnitude larger 
than the DAMA Nal(Tl) background in the same energy region) by hardware vetoes, 

k In fact, no direct measurement performed with neutron sources or generators has been reported 
up to now, although several bolometers have been irradiated with neutrons along the past decade. 
For the sake of completeness, we remind that a measurement of the response of a TeC>2 bolometer 
to surface 224 Ra recoiling nuclei has been reported in ref. EE1; this measurement, although its 
importance, does not represent a determination of the quenching factor for a bolometer since neither 
the target-nuclei nor the whole bulk of the detector are involved. Moreover, these values cannot of 
course be extended to whatever kind of bolometer. 

'For comparison we remind that the data of the DAMA/LXe set-up, which has a similar light 
response, are analysed by using the much more realistic and safer software energy threshold of 13 
keV^. 



Table 1: Features of the DAMA/Nal results on the WIMP annual modulation signature during the 
first four annual cycles (57986 kg x day exposure) 20 > 22 > 23 > 24 > 25 > 26 > 27 > 28 ) anc j those of refs. 39 > 40 > 41 ). 





DAMA/Nal 


CDMS-I 


Edelweiss-I 


Zeplin-I 


Signature 


annual 
modulation 


None 


None 


None 


Target-nuclei 


23 Na, 127 I 


natQ e 


"°'Ge 


nct Xe 


Technique 


well known 


poorly 
experienced 


poorly 
experienced 


critical optical 
liquid/gas interface 
in this realization 


Target mass 


~ 100 kg 


0.5 kg 


0.32 kg 


~3kg 


Exposure 


57986 kg X day 


15.8 kg X day 


8.2 kg X day 


280 kg X day 


Depth of the 
experimental site 


1400 m 


10 m 


1700 m 


1100 m 


Software energy 
threshold 
(electron 
equivalent) 


2 keV 
(5.5 - 7.5 p.e./kcV) 


10 keV 


20 keV 


2 keV 
(but: cr/E = 100% 
mostly 
1 p.e./keV;B3) 

(2.5 p.e./keV 
for 16 days; El) 


Quenching 
factor 


Measured 


Assumed = 1 


Assumed = 1 


Measured 


Measured event 
rate in low 
energy range 


~ 1 cpd/kg/keV 


~ 60 cpd/kg/keV 
(10 5 events) 


2500 events 
total 


~ 100 cpd/kg/keV 


Claimed events 
after rejection 
procedures 




23 in Ge, 4 in Si, 
4 multiple evts in Ge 
+ MonteCarlo on 
neutron flux 





~ 20-50 cpd/kg/keV 

after rejection and 
?? after standard PSD 
1411491 


Events satisfying 
the signature 
in DAMA/Nal 


modulation 
amplitude 
integrated over the 
given exposure 
~ 2000 events 








Expected number 

of events from 
DAMA/Nal effect 




from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


depends on 
the models 
(even zero) 



by fiducial volume cuts and, largely, by applying down to few keV a standard pulse 
shape discrimination procedure, although the LXe scintillation pulse profiles (pulse 
decay time < 30 ns) are quite similar even to noise events in the lower energy bins 
and in spite of the poor light response. Quantitative informati on on experimental 
quantities related to the used procedures has not yet been given EES! 



In addition to the experimental aspects, these experiments - which cannot per- 
form any model independent comparison with the result of DAMA/Nal - generally 
uncorrectly/partially quote the results of published quests for a purely SI coupled 
candidate in given model frameworks and ignore the effects of the relevant uncertain- 
ties in the astrophysical, nuclear and particle physics assumptions and experimen- 
tal/theoretical parameters values taken by each experiments as well as the interpre- 
tation of the DAMA/Nal model independent effect in terms of candidates with other 
kind of coupling. All that can significantly vary the result of any comparison (even 
when assuming as correct the evaluation of the selected nu mber of even ts, the energy 
scale and the energy threshold determinations given in refs. EHEIIEEIIl an( } fa e US ed 
efficiencies, etc.). In addition, there exist scenarios to which Na and I are sensitive 
and other nuclei, such as e.g. nat Ge, na *Si and nat Xe, are not. 

In conclusion: 

1. no other experiment, whose result can be directly comparable in a model inde- 
pendent way with that of DAMA/Nal, is available so far. 

2. as regards in particular CDMS-I, EDELWEISS-I and Zeplin-I, e.g.: 

i) they are insensitive to the model independent WIMP annual modulation 
signature exploited by DAMA/Nal; ii) they use different methodological ap- 
proaches, which do not allow any model independent comparison and they have 
different sensitivities to WIMPs; in particular, the number of counts they could 
expect on the basis of the model independent DAMA/Nal result varies from few 
to zero events depending on the models, on the assumptions and on the the- 
oretical/experimental parameters' values adopted in the calculations; iii) they 
make neither correct nor complete comparisons with the DAMA/Nal experi- 
mental result; iv) they use extremely poor statistics; v) they reduce their huge 
measured counting rate of orders of magnitude by various rejection procedures 
claiming for very optimistic rejection powers; vi) their energy scale determina- 
tion and/or energy threshold appear questionable (in the first two cases because 
of the quenching factors values and in the second because of the poor number 
of photoelectrons/keV); etc. 

5.1.3.2 ... from indirect searches 

As regards results from indirect searches, model dependent analyses investigat- 
ing possible up-going muons arising from WIMP annihilation in celestial bodies have 
been carried out by large experiments deep underground such as e.g. Macro and Su- 
per kamiokande. The comparison is strongly model dependent; anyhow it has already 
been shown in ref. that even in the simple model frameworks considered in ref. 

the model dependent limit on up-going muon flux obtained by Macro (and thus 
that by Super kamiokande, which is only slightly more stringent) could cut only part 
of the allowed configurations in MSSM. 



On the other hand, as regards antimatter searches carried out outside the atmo- 
sphere, we remind the analysis of the HEAT balloon-borne experiment performed in 
ref. where an excess of positrons with energy ~ 5 — 20 GeV has been found and 
has been interpreted in terms of WIMP annihilation. We also mention the analysis 
of ref. which already su gges ts the presence of a 7 excess from the center of the 
Galaxy in the EGRET data^* which matches with a possible WIMP annihilation 
in the galactic halo and is not in conflict with the DAMA/Nal model independent 
result. 

We take this occasion to stress, however, that the specific parameters of a WIMP 
candidate (mass and cross sections), which can be derived from the indirect searches, 
critically depend on several assumptions used in the calculations such as the esti- 
mation of the background, the halo model, the amount of WIMP in the galactic 
dark halo, the annihilation channels, the transport of charged particle to the Earth, 
etc.; thus, they have the same relative meaning as those obtained in the quest for a 
candidate in direct search approach as previously described. 

5.1.3.3 Conclusions 

In conclusion, no model independent comparison with the DAMA/Nal effect is 
available. Only few model dependent approaches have been used in direct searches to 
claim for a particular model dependent comparison, which appears in addition - as 
mentioned above - neither based on solid procedures nor fully correct nor complete. 
On the other hand, the indirect search approaches, which can also offer only model 
dependent comparisons, are not in contradiction or in substantial agreement with the 
DAMA/Nal observed effect. 

Thus, the interest in the further available DAMA/Nal data is increased and the 
analysis of the cumulative data of the seven annual cycles will be released in near 
future. 

6. The new DAMA /LIBRA set-up 

After the completion of the data taking of the DAMA/Nal set-up, the procedures 
to install the new LIBRA set-up have been carried out. In particular, improvements 
have been realized in the experimental site, in the Cu box, in the shield and in the 
available detectors. 

The LIBRA set-up consisting of ~ 250 kg of radiopure Nal(Tl) is made of 25 
detectors, 9.70 kg each one. The new detectors have been realized thanks to a second 
generation R&D with Crismatec company, by exploiting in particular new radiop- 
urification techniques of the Nal and Til powders. In the framework of this R&D 
new materials have also been selected, prototypes have been built and tested and a 
devoted protocol has been fixed. 

The full dismounting of the ~ 100 kg Nal(Tl) set-up, the improvements mentioned 



above and the installation of all the detectors including a new PMTs' shield have 
been completed at end 2002. All the related procedures have been performed in 
HP Nitrogen atmosphere by using special masks connected to air bottles to avoid 
that the inner part of the Cu box, the detectors, the new PMTs' shield, etc. would 
be in contact with environmental air, that is to reduce at most possible surfaces' 
contamination by environmental Radon. Some pictures taken during the installation 
of the LIBRA detectors are shown in Fig. ^2 Before this installation, all the Cu 
parts have been chemically etched following a devoted protocol and maintained in 
HP Nitrogen atmosphere until the installation. 




Figure 11: Left picture: during the detectors installation in HP Nitrogen atmosphere. Right picture: 
view at end of the detectors installation. All the used materials have been deeply selected for 
radiopurity (see for example the cables with teflon envelop). 

LIBRA will offer a relevant competitiveness e.g. because of its: i) high sensitiv- 
ity; ii) standard and well defined operating procedures; iii) well known technology; 
iv) proved possibility of an effective control of the experimental conditions during 
several years of running; v) high duty cycle; vi) possibility to deeply investigate the 
WIMP model independent signature; vii) sensitivity to both SI and SD couplings; 

viii) favoured sensitivity in some of the possible particle and astrophysical models; 

ix) high benefits/cost. 

The main aim of this new LIBRA set-up is to further investigate the WIMP 
component in the galactic halo with increased sensitivity thanks to the larger exposed 
mass and to the higher overall radiopurity. Moreover, when applicable, it offers 
also pulse shape discrimination capability and the possibility to achieve competitive 
results in the searches for other rare processes. For example, it would reach in a 
relatively short time a sensitivity of ~ 3 x 10 27 y for possible Pauli-exclusion-principle 
violating processes, a sensitivity of ~ 10 24 — 10 25 y for possible charge-non-conserving 
processes in 23 Na and in 127 I (depending on the counting rate in the energy region of 
interest), a sensitivity of ~ 10 27 y for the nucleon and di-nucleon decay into invisible 
channels, a sensitivity of ~ 10~ 10 GeV -1 on the axion-photon coupling constant from 
the investigation of solar axions; in addition, it could explore mass for exotic Dark 



Matter candidates, such as e.g. the SIMPs, up to ~ 10 17 GeV and neutral nuclearites 
flux of ~ 5 x 10~ 12 s _1 cm _2 sr _1 . 

7. Conclusion 

In this paper recent results achieved by the DAMA experiment at the Gran Sasso 
National Laboratory of I. N.F.N, have been summarized. In particular, DAMA/Nal 
has been a pioneer experiment running at LNGS for several years and investigating 
as first the WIMP annual modulation signature with suitable sensitivity and control 
of the running parameters. It has pointed out the presence of a modulation satisfying 
the several peculiarities of a WIMP induced effect and the absence of any possible 
systematic effects or side reactions able to mimic it, reaching a significant model 
independent evidence. As a corollary result, it has also pointed out the complexity of 
the quest for a WIMP candidate because e.g. of the present poor knowledge on the 
many astrophysical, nuclear and particle physics aspects and on its nature. 

After the completion of the data taking of the ~ 100 kg Nal(Tl) set-up (on 
July 2002; full statistics of 107731 kg • d in progress to be released), as a result 
of the continuous efforts toward the creation of ultimate radiopure set-ups, the new 
DAMA/LIBRA has been realized and installed. Further R&D efforts for ultimate 
radiopurification of Nal(Tl) detectors are also starting again. 
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